Dimensional tuning of electronic states under strong and frustrated interactions 
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We study a model of strongly interacting spinless fermions on an anisotropic triangular lattice. At half- 
filling and the limit of strong repulsive nearest-neighbor interactions, the fermions align in stripes and form an 
insulating state. When a particle is doped, it either follows a one-dimensional free motion along the stripes or 
fractionalizes perpendicular to the stripes. The two propagations yield a dimensional tuning of the electronic 
state. We study the stability of this phase and derive an effective model to describe the low-energy excitations. 
Spectral functions are presented which can be used to experimentally detect signatures of the charge excitations. 
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Electronic properties of metals are essentially understood 
on the basis of band theory which straightforwardly reflects 
the actual crystal geometry. In certain strongly correlated ma- 
terials, however, the electronic interactions can significantly 
modify the spacial geometry of the wave functions (metallic- 
ity). For example, in high-T c cuprates, the charges on the 
square lattice may align in stripes and the metallic proper- 
ties become one dimensional (ID) QJJJ . Also in manganites, 
the strong correlations of spin and orbital degrees of freedom 
lead to orbital orders in which the metallicity is confined to 
a particular spatial direction [2]. Here the density of states 
reflects correlation effects and causes an effective modulation 
of the lattice geometry and thus leads to a change of the in- 
trinsic thermodynamical and transport properties. A recently 
proposed way to tune the effective geometry of a lattice is to 
introduce frustrations in the electronic interactions. This is 
seen in a model of spinless fermions (i.e., fully spin polarized 
electrons) on a triangular lattice with strong nearest neighbor 
repulsion which is referred as a "pinball liquid". Part of the 
fermions form a Wigner crystal on one of the three sublattices. 
The other fermions form a metallic state on the remaining sites 
along a honeycomb lattice (3D. 

In this Letter, we investigate a tuning of the effective di- 
mensionality of the metallicity due to partially frustrated and 
strong interactions. We focus on an anisotropic triangular 
lattice system and start from an insulator at half-filling with 
striped charge order. Doping of particles yields a ID metal- 
lic state as in the case of the square lattices in cuprates. In 
addition, we find another non-trivial propagation in which the 
doped charges fractionalize and change the stripe geometry. 

We focus only on the charge degrees of freedom and start 
from a basic t-V model Hamiltonian which reads 
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Here cj (c-) are creation (annihilation) operators of fermions 

and rij (=c]cj) are number operators. The interactions act 
only between neighboring sites (ij). Anisotropies of the 
hopping amplitudes and repulsion strengths are given by 
(Uj,Vij) = (t',V) > (0, 0) for vertical bonds (along the 



2/-axis in, e.g., Fig. [U and (t, V) > (0, 0) for the remaining 
bond directions. We start from the limit of strong anisotropic 
interactions, i.e., t/V = 0, t/V = 0, t/\V - V\ = 0, and 
assume for simplicity that t' = t > 0. The ground state at 
half-filling has an insulating charge order which minimizes 
the repulsive nearest-neighbor interactions [4] . We consider 
the effect of doping either a particle or a hole into the insulat- 
ing state. Depending on the anisotropy of the interactions, 
we find two different scenarios. The case V < V is de- 
picted schematically in Fig. Uta). The ground state at half- 
filling has a vertical stripe-order. An added particle can move 
freely along a line of empty sites without changing the po- 
tential energy of 4V. All processes of moving orthogonal to 
the stripes are prevented by an increase of potential energy by 
2(V — V). This implies that the added particle is confined to 
a ID motion along the ?/-axis with a dispersion —2tcos(k y ) 
for wave-number k y . Doping of one hole decreases the po- 
tential energy by 2V' and yields a spectrum identical to that 
of the added particle. For the case V > V, the ground state 
at half filling has a high degeneracy. The ground states are 
formed by isolated configurations which have stripes perpen- 
dicular to the anisotropic bonds with a 2 Nx -fold degeneracy 
Q4Q . Figure Hlb) shows the horizontal stripe configuration as 
one example of degenerate ground state configurations. An 
added particle can propagate along the stripes as in the for- 
mer case. However, an additional propagation can take place 
now. The doped particle interacts with the stripes in the verti- 
cal direction by two V -bonds. These two bonds can separate 
orthogonally to the stripes without increasing the potential en- 
ergy. Each of them is carrying a charge of — e/2. These local 
defects are in the following regarded as fractional charges. As 
long as we remain in the limit of strong interactions, i.e., as- 
sume t/V = 0, the fractional charges are deconfined. When a 
particle is removed, a similar picture arises from the fraction- 
alization of a hole. 

In the above discussion, we have neglected any quantum 
fluctuation in t/V. Next, we will investigate how quantum 
fluctuations affect the ground state and calculate the low- 
lying excitations. We numerically diagonalize the Hamilto- 
nian in Eq. (Q]) exactly on finite size clusters and derive a 
strong-coupling effective model by a perturbative approach. 
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We start by calculating the local density of states (DOS), 
D± (v) = Y^k A ± (k,u) of the full Hamiltonian in Eq. £Q) 
by means of the Lanczos continued fraction method [5]. Here, 
A ± (k, uj) is the spectral function from adding (+) and remov- 
ing (— ) a particle. We use aN — N x x N y = 4x6 site cluster 
with periodic boundary conditions. Figs. [2a) and[2b) show 
the DOS for V < V and V > V in the strong coupling re- 
gion with V + V = 2Vq = 300. The spectrum of the first 
case (V < V) in Fig. 0a) has a width of At and corresponds 
exactly to a ID dispersion in a single-particle picture. For 
the latter case (V > V), we find instead an incoherent broad 
structure in Fig. 0(b). The bandwidth is «8£ and thus twice 
as big as in the former case. This indicates the existence of a 
collective excitation of fractionalized charges. 

In the following, we focus on the interesting case in which 
fractionalization of charge occurs. We begin by estimating the 
parameter range for which the fractional-charge phase sus- 
tains. In this phase, the ground-state energy depends on the 
interaction in a specific manner, i.e., we find E ~ — 
(2 V + 4t) and \N + (2V + 2 V - 4t) for one hole and one 
particle doping, respectively. The first term originates from 
the stripes and the second term is the interaction energy of the 
doped charges. On the other hand, a pinball liquid phase is 
found for V' ~ V (i.e., Vq ~ V) which also shows linear In- 
dependence in energy as discussed in Ref . [4] . We distinguish 
the phases by their different linear slopes in the energy E(V). 
The phase boundaries are then identified for up to two parti- 
cle/hole doping (~ 8% doping) as shown in Fig. 0c). From 
these results we can expect that the fractional-charge phase re- 
mains stable for a large range of anisotropic interactions even 
if the doping level is increased. Note that if we increase the 
system size towards the bulk limit, the phase boundary for the 
doped particle number of 0(1) will asymptotically approach 
the solely ^-dependent line at V — Vo = — ^ • We thus con- 
clude that the physics we discuss is relevant for a large range 
of interaction strength and different doping levels. 

Next, in order to understand in detail how the fluctua- 
tions affect the low-energy excitation, we perform a pertur- 
bative treatment in the limit of strong couplings. We be- 
gin with the half-filled case which has a degenerate ground- 
state manifold of stripe configuration and include the pertur- 
bative processes in t/V, t/V and t/\V — V\ up to fourth 
order. The first contribution which lifts the degeneracy are the 
fourth-order exchange processes around plaquettes as shown 
in Fig.Oa). We distinguish between (A) plaquettes with odd 
number of particles (i.e., either one or three) and (B) pla- 
quettes with two particles. These processes give energy cor- 
rections of —64 and +§64 per plaquette, respectively, where 
e 4 = 2t A /((2V - V) 2 V) > 0. The different signs originate 
from an exchange of fermions. The conversion of one plaque- 
tte from (A) to (B) increases the energy by A = |e4. Then, 
at half-filling, the afore mentioned degeneracy of stripes at 
t/V = is lifted and the horizontal stripe filled with (A) pla- 
quettes becomes a unique ground state. The lowest excitation 
energy is N y A which results from a collective rearrangement 
of charges to another stripe |4j]. In the fractional-charge phase 



(a) vertical stripe+1 electron V'<V 
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FIG. 1: (Color online) Representative ground-state configurations of 
the Hamiltonian in Eq.(l) on the anisotropic triangular lattice in the 
limit of strong correlations (t/V = 0) with one particle doped at 
half-filling. Panel (a) show the case on a vertical stripe order at 
V' < V where an added particle can move only along the stripes 
of empty sites. Panel (b) shows the horizontal stripe at V' > V, 
where the fractionalization is illustrated by the separation of the two 
bold vertical bonds in the ^/-direction. 
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FIG. 2: (Color online) Numerical results on a N x x N y = 4x6 = 
24-site cluster. DOS of Eq. (1) with t = t' = 1 for (a) (V, V) = 
(200, 100) and (b) (V, V) = (100, 200). Panel (b) compares the 
DOS of Eq. (1) for V' = 200, V = 100 (bold line) with the effective 
model (shaded region) of Eq. (2) on top of each other, (c) shows 
the ground state phase diagram. The phase boundary is shown for 
different doping levels. In the bulk system the phase boundary at low 
doping limit approaches the broken line, where we have fractional- 
charge phase at V — Vo < 
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near half-filling, those plaquette-exchange processes lead to 
a confinement. This can be seen directly in Fig. Ob); the 
separation of fractionalized charges converts (A) plaquettes 
to (B) plaquettes. The confinement potential A c (|yi — 2/2 1)> 
which sums up all the possible perturbative processes up to 
fourth order (excluding a constant shift), is shown in Fig. [3c). 
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(A)-MB) A =+ -| £ 4 /plaquette 



FIG. 3: (Color online) (a) Fourth order ring exchange pro- 
cesses along the plaquettes; (A)/(B) lowers/raises the energy by 
(— €4)/(+f 64) (both the clockwise and anti-clockwise processes are 
counted), (b) The case when one particle is added. The shaded 
plaquettes between the two vertical bonds are the plaquettes re- 
placed from (A) to (B) when the bonds are separated by the distance 
I yi — V2 1 . (c) Energy correction A c including all the perturbation 
processes up to fourth order as a function of | y± — yi | for several 
parameter choices. The gradient is 2 A = 5 £4, which comes from 
the conversion of (A) to (B) plaquettes. 

Here y\ and 1/2 are the ^/-components of the locations of the 
two fractions. The main feature is a linear contribution of 
4 A per two lattice spacings (for 1 2/1 — 2/2 1 > 3) which arises 
from part of the fourth-order processes where the conversion 
of (A) to (B) plaquettes takes place. The above findings al- 
low us to introduce the following effective Hamiltonian Il4ll 
to describe the low energy excitations of the doped system for 

^Pr = E V (Vh" + HX ') V + A ( N W - N W)' < 2 > 

(i,3) 

The projector V projects onto the manifold of allowed con- 
figurations which minimize the interaction term of Eq. (1). 
The projected configurations are connected by first order in 
Uj . The second term takes into account the confinement ef- 
fects which result from the conversions of (A) to (B) plaque- 
ttes. The operators N^m and count the number of cor- 
responding plaquettes. We validate the effective Hamiltonian 
by comparing the obtained electron doped DOS with the one 
obtained from the Full Hamiltonian Eq. (Q]). As presented in 
Fig. Etb), the DOS of Eq. © for (V, V) = (100,200) in 
solid line agrees perfectly well with the one from the effec- 
tive model on the same cluster presented as a shaded region. 
The fact that we have to consider only the manifold of stripe- 
based configurations reduces the size of the Hilbert space of 
TLp r drastically as compared to the full Hilbert space. Thus 
we can perform a careful numerical analysis of the effects of 
the confinement potential. 

Let us first focus on the case, A « 0. When a particle is 
added, it moves either along the stripes freely by t or frac- 
tionalizes vertically to the stripes by t' . This is observed in 
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FIG. 4: (Color online) Panel (a) shows the energy dispersion of 
the effective Hamiltonian Eq. (0 for a small confining potential 
4A = 10~ 6 on a N = 8 x 12 cluster. Panel (b) and (c) show the 
local density of states D(oo) for different confining potentials and for 
different anisotropics in the hopping amplitudes, respectively. We 
use for the spectrums the N = 4 x 32 or 8 x 24 cluster which are 
large enough to give qualitatively the same results (we checked the 
finite size effect up to N y = 42). Arrows in panel (b) indicate the 
three intrinsic peak positions. The energy is shifted by a constant 
C(A) for a clear presentation. Inset: Schematic illustration of the 
variation of the spectrums. Arrows denote the (non-fractionalized) 
van Hove singularity. 



the dispersions within the finite cluster. Figure [4]^) shows 
the two-dimensional dispersion for t = t' = 1. The disper- 
sions along r — X are nearly flat since the x-direction mo- 
tion is sacrificed by the increase of coherence of fractional 
charges in ^-direction. There are dispersive degenerate modes 
with a width of 4t along M-Y. These modes originates from 
the non-fractionalized portion which can only propagate in in- 
direction. They will be shown later to contribute to the weak- 
ened van Hove singularity. The feature of the fractionalization 
in the ^/-direction is detected along the X-M line. We find 
densely packed branches which are following a cosine disper- 
sion. The bandwidth is w 8t' and is almost independent of N y 
while the number of modes is N y /2. Thus, the dispersions ex- 
trapolate towards a continuum in the bulk limit. The existence 
of continuum has already been discussed in terms of domain 
wall excitations in ID Wigner lattices and is visible in the ex- 
act solution 0]. The two fractions share the total momentum 
k with arbitrary ratio so that at each /c -point the energy can 
take almost any values of energy within a certain range. 

The effect of the confinement potential becomes clear by 
comparing the DOS obtained for different choices of A (see 
Fig. EJb)). If the confinement is small, we find a three peak 
structure which results from the superposition of the ID van 
hove singularities at uo = ±2t and a uo = 0-centered broad 
single peak structure. The former singularity is a result of 
the non-fractionalized propagation in the ^-direction. The lat- 
ter single peak is due to the fractionalization of charge and 



4 



is consistent with the previous ID study |6]. The weight of 
the single peak is suppressed by a very large confining poten- 
tial A. However, even for confinement potentials as large as 
A ~ t, we still find broad incoherent structure expanding 
over the energy range of > 2t. This is a clear indica- 
tion that although the two fractions cannot separate to infinite 
distance under confinement, the coherence due to t' is still 
large and they form bound pairs with large spatial extent (usu- 
ally of order t'/A in units of the lattice constant). A simi- 
lar transformation of the DOS is observed when we consider 
the case without confinement potential but change instead the 
anisotropic hopping parameters by hand. Figure HJc) shows 
the electron doped spectrum for different ratios of t' jt. For 
t = 0, t' 7^ 0, the fractional charges form a broad single 
peak structure. When introducing a finite t, the two ID van- 
Hove singularity peaks appear, and it finally becomes a pure 
ID DOS of free fermion along the ^-direction at t' = 0, t ^ 0. 

In summary, we find a state with exotic quasiparticles de- 
scribed by a t- V model on the triangular lattice in the presence 
of strong and anisotropic interactions (V > V ^> t, t'). An 
added particle or hole decays into bound pairs of fractional 
charges with large spatial extent. Not only the added particle 
but also the other particles move cooperatively and show col- 
lective excitation behavior in the anisotropic direction. Even 
though the system is originally two dimensional, the strong 
and anisotropic interactions first lead to a ID like charge or- 
dering at half-filling. The dispersion of the quasiparticle is the 
combination of the ID free (t) and ID collective it') propa- 
gations along and perpendicular to stripes, respectively. The 
effective dimension is tuned from one to two-dimension de- 
pending on the anisotropy oft'/t, on the perturbative interac- 
tions and on the original interaction strength. It is remarkable 
that such delicate and non-trivial dimensional tuning occurs in 
such a simple model. 

Fermions and bosons on certain frustrated lattices have 
shown to behave in a similar manner, e.g., a pinball liquid 
found previously in the present model L3] and a supersolidity 
in a hard core bosonic model [7] can be compared. We find 
that in the present system, a similar hard-core bosonic model 
has a diagonal stripe ground state at half-filling instead of the 
horizontal stripe and it shows similar fractionally charged ex- 
citations. These fractionalized excitations are bound by a con- 
finement potential of one-half the magnitude of that of the 
fermionic model. The different details originate simply from 
the different statistics, i.e., the absence of a fermionic sign 
when exchanging two hard-core bosons. 

Finally we note that such dimensional tuning is not limited 
to the anisotropic triangular lattice. We find similar effects 
also on the anisotropic kagome lattice at around 2/3 -filling. 
Such picture in between the classical order and the quantum 
liquid might more generally be expanded to several other frus- 
trated lattices as well. 

The two-dimensional charge ordered stripe state originating 
from strong Coulomb interactions is already an established 
phenomenon in organic solids 0]. To search for its further 



interesting possibilities, optics become a powerful tool which 
aim for the photo-induced phase transitions as in EDO-TTF 
Jl], 6>-ET 2 X. |H and a-ET 2 I 3 LH]. We expect the highly 
conducting photo-induced state based on the horizontal stripe 
in the latter two materials to be a good candidate to realize 
the scenario presented in this letter. Fractional charge itself 



has been discussed in many theoretical articles I& LLJ, 11311 . 
Although the present quantum fractionalization due to frus- 
tration is not exactly the one found in the classical particle 
picture, it is certainly exotic and is the strong indications for a 
collective many body effect. Its experimental proof should be 
awaited. 
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